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Gradient sensing requires at least two measurements at different points in space. These measure¬ 
ments must then be communicated to a common location to be compared, which is unavoidably 
noisy. While much is known about the limits of measurement precision by cells, the limits placed 
by the communication are not understood. Motivated by recent experiments, we derive the funda¬ 
mental limits to the precision of gradient sensing in a multicellular system, accounting for commu¬ 
nication and temporal integration. The gradient is estimated by comparing a “local” and a “global” 
molecular reporter of the external concentration, where the global reporter is exchanged between 
neighboring cells. Using the fluctuation-dissipation framework, we find, in contrast to the case when 
communication is ignored, that precision saturates with the number of cells independently of the 
measurement time duration, since communication establishes a maximum lengthscale over which 
sensory information can be reliably conveyed. Surprisingly, we also find that precision is improved 
if the local reporter is exchanged between cells as well, albeit more slowly than the global reporter. 

The reason is that while exchange of the local reporter weakens the comparison, it decreases the 
measurement noise. We term such a model “regional excitation-global inhibition” (REGI). Our 
results demonstrate that fundamental sensing limits are necessarily sharpened when the need to 
communicate information is taken into account. 

Cells sense spatial gradients in environmental chemi¬ 
cals with remarkable precision. A single amoeba, for ex¬ 
ample, can respond to a difference of roughly ten attrac- 
tant molecules between the front and back of the cell [1] . 

Cells are even more sensitive when they are in a group: 
cultures of many neurons respond to chemical gradients 
equivalent to a difference of only one molecule across an 
individual neuron’s axonal growth cone [2]; clusters of 
malignant lymphocytes have a wider chemotactic sensi¬ 
tivity than single cells [3]; and groups of communicat¬ 
ing epithelial cells detect gradients that are too weak 
for a single cell to detect [4]. More generally, collec¬ 
tive chemosensing properties are often very distinct from 
those in individual cells IlIMZ]- These observations have 
generated a renewed interest in the question of what sets 
the fundamental limit to the precision of gradient sens¬ 
ing in large, spatially extended, often collective sensory 
systems. 

Fundamentally, sensing a stationary gradient requires 
at least two measurements to be made at different points 
in space. The precision of these two or more individ¬ 
ual measurements bounds the gradient sensing precision 
In their turn, each individual measurement is lim¬ 
ited by the finite number of molecules within the detector 
volume and the ability of the detector to integrate over 
time, a point first made by Berg and Purcell (BP) [5]. 

More detailed calculations of gradient sensing by spe¬ 
cific geometries of receptors have since confirmed that 
the precision of gradient sensing remains limited by an 
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expression of the BP type [TOHI^ . 

However, absent in this description is the fundamental 
recognition that in order for the gradient to be measured, 
information about multiple spatially separated measure¬ 
ments must be communicated to a common location. 
This point is particularly evident in the case of multi¬ 
cellular sensing: if two cells at either edge of a popu¬ 
lation measure concentrations that are different, neither 
cell “knows” this fact until the information is shared. 
This is also important for a single cell: information from 
receptors on either side of a cell must be transported, 
e.g., via diffusive messenger molecules, to the location of 
the molecular machinery that initiates the phenotypic re¬ 
sponse. How is the precision of gradient sensing affected 
by this fundamental communication requirement? 

As discussed recently in the context of instantaneous 
measurements [1] , the communication imposes important 
limitations. First, detection of an internal diffusive mes¬ 
senger by cellular machinery introduces its own BP-type 
limit on gradient sensing. Since the volume of an in¬ 
ternal detector must be smaller than that of the whole 
system, and diffusion in the cytoplasm is often slow, such 
an intrinsic BP limitation could be dramatic. Second, in 
addition to the detection noise, the strength of the com¬ 
munication itself may be hampered over long distances 
by messenger turnover. This imposes a finite lengthscale 
over which communication is reliable, with respect to the 
molecular noise. However, the communicating cells can 
integrate the signals over time [5], improving detection 
of even very weak messages. Whether such integration 
can lift the communication constraints has not yet been 
addressed. 

To analyze constraints on gradient sensing in spatially 
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extended systems with temporal integration, we use a 
minimal model of collective sensing based on the lo¬ 
cal excitation-global inhibition (LEGI) approach [T3]. 
This sensory mechanism uses a local and a global inter¬ 
nal reporter of the external concentration, where only 
the global reporter is exchanged and averaged among 
neighboring cells. Comparison of the two reporters then 
measures if the local concentration is above or below 
the average, and hence whether the cell is on the high- 
concentration edge of the population. We analyze the 
model using a fluctuation-dissipation framework m , to 
derive the precision with which a chemical gradient can 
be estimated over long observation times. In the case 
where the need to communicate is ignored, the preci¬ 
sion would grow indefinitely with the number of cells. In 
contrast, we find that, communication imposes limits on 
sensing even for long measurement times. Furthermore, 
the analysis reveals a counterintuitive strategy for opti¬ 
mizing the precision. We find that if the local reporter is 
also exchanged, at a fraction of the rate of the global re¬ 
porter, the precision can be significantly enhanced. Even 
though such exchange makes the two compared concen¬ 
trations more similar, which weakens the comparison, it 
reduces the measurement noise of the local reporter. This 
tradeoff leads to an optimal ratio of exchange rates that 
maximizes sensory precision. We discuss how the predic¬ 
tions of our analysis could be tested experimentally. 


I. RESULTS 

We discretize the spatially-extended gradient sensor 
into compartments. These can be whole cells or their 
parts, but we will refer to them as cells from now on. 
There is a diffusible chemical whose concentration varies 
linearly in space. The chemical gradient defines a direc¬ 
tion within the sensor, and we focus on a chain of cells 
along this direction (Fig. Numbering the cells from 
n = 1 to each cell experiences a local concentration 
Cn = cn — (N — n)ag, where cn is the background con¬ 
centration, g is the concentration gradient, and a is the 
linear size of each cell. We choose without loss of gen¬ 
erality to have g > 0 and to reference the background 
concentration at cell N, which is then at the higher edge 
of the gradient. We focus on this cell because we imagine 
it will be the first to initiate a phenotypic response, such 
as proliferating or directed motility. Finally, we focus on 
the limit ag/cN 1 because limits on the sensory preci¬ 
sion will be the most important for such small, hard to 
measure gradients. 


A. An idealized detector 

First, we consider the case when the two edge cells 
form an idealized detector, in the sense that each cell 
counts every external molecule in its vicinity, and one 
cell “knows” instantly and perfectly the count of the 
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FIG. 1: Spatially extended gradient sensing. (A) A 

chain of N compartments or cells is exposed to a linear pro¬ 
file of a diffusible chemical. (B) In an idealized detector, the 
two edge cells communicate their measurements perfectly and 
instantly. (C) In our model, bound receptors (R) activate a lo¬ 
cal (X) and global reporter molecules (Y), and Y is exchanged 
between cells for the communication. 


other (Fig.[^). The gradient could then be estimated by 
the difference in the concentration measurements made 
by the two cells IS M- The mean of this difference 
is A = Cat — Cl = {N — l)ag, and its error is given 
by the corresponding errors in the two measurements, 
(i5A)^ = ((5ci)^ -I- (Scn)"^, under the assumption that the 
measurements are independent. 

Berg and Purcell showed [S] that the fractional er¬ 
ror in each measurement is not smaller than (Sc/c)^ ^ 
l/{aDTc), where D is the diffusion constant of the lig¬ 
and, and T is the time over which the measurement is 
integrated. This expression has an intuitive interpreta¬ 
tion: the fractional error is at least as large as the Poisson 
counting noise, which scales inversely with the number 
of molecular counts. The number of counts that can be 
made in a time T is given by the number of molecules in 
the vicinity of a cell at a given time, roughly ca^, multi¬ 
plied by the number of times the diffusion renews these 
molecules, T/r, where t ^ /D. This product is aDTc. 

The error in the gradient estimate is then given by 
(i5A)^ ~ ci/(aDT) + cjq/{a-DT). For sufficiently small 
gradients, such that ci « cat, this becomes 


Thus, for an idealized detector, the error in the gradient 
estimate is limited entirely by the error in the measure¬ 
ments made by each of the edge cells. In general, we 
could think of the measurement at either edge being per¬ 
formed by a region that is larger than a single cell. Since 
each region could not be larger than the whole system, 
the highest precision is obtained when ~Aa replaces a 
in Eqn. [l] This result has been derived more rigorously 
m, and apart from a constant prefactor, Eqn. [^indeed 
provides the estimation error in the limit of large detector 
separation and fast detection kinetics. More complex ge¬ 
ometries, such as rings of detectors [10], or detectors dis¬ 
tributed over the surface of a circle m or a sphere 
have also been considered, and Eqn. again emerges as 
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the corresponding bound, with the lengthscale a replaced 
by that dictated by the specific geometry. 

Eqn. can be combined with the mean A to produce 
the signal-to-noise ratio (SNR) for gradient detection 

1 __ ( 2 ) 

SNR \aJ [{N- l)ag]'^aDT '' 

This expression again has a clear interpretation: the SNR 
increases if the external molecules diffuse more quickly 
{D) or are more sharply graded (g), or if the detectors are 
larger (a), are better separated {N), or integrate longer 
(T). Yet, the SNR is worse for a larger background con¬ 
centration (ctv), since it is more difficult to detect a small 
gradient on top of a larger background. The measures de¬ 
fined in Eqns.j^andj^are conceptually equivalent only in 
the case of low background concentration, when the dif¬ 
ference A is comparable to the background concentration 
Cat (i. e. Ci « 0). While much of the field has focused on 
Eqn. [2 here we are concerned with the opposite case: the 
fundamental limits to the detection of small gradients on 
a large background. Therefore from here on we focus on 
the SNR, Eqn. § 

B. Accounting for the need to communicate 

Eqn. cannot be a fundamental limit because it ne¬ 
glects a critical aspect of gradient sensing: the need to 
communicate information from multiple detectors to a 
common location. Indeed, the idealized detector implies 
the existence of a “spooky action at a distance” [TS], i. e. 
an unknown, instantaneous, and error-free communica¬ 
tion mechanism. What are the limits to gradient sensing 
when communication is properly accounted for? 

To answer this, a model of gradient sensing must be 
assumed. A naive model would allow each cell access to 
information about the input measured and broadcast by 
every other cell. This would require a number of private 
communication channels that grows with the number of 
cells, which is not plausible. A realistic alternative that 
would involve just one message being communicated is for 
each cell to have access to some aggregate, average infor¬ 
mation, to which all comparisons are made. There are 
a few such models [TBHIH], and our choice among them 
is guided by the fact that collective detection of weak 
gradients is observed in steady state, and over a wide 
range of background concentration in both neurons [ 2 ] 
and epithelial cells [4] . This supports an adaptive spatial 
(rather than temporal) sensing, such as can be imple¬ 
mented by the local excitation-global inhibition (LEGI) 
mechanism m- 

The LEGI model is illustrated in Eig. [Tp. Each 
cell contains receptors that bind and unbind external 
molecules with rates a and /i, respectively. Bound re¬ 
ceptors (R) activate both a local (X) and a global (Y) 
intracellular species with rate /3. Deactivation of X and 
Y occurs spontaneously with rate v. Whereas X is con¬ 
fined to each cell, Y is exchanged between neighboring 


cells with rate 7 , which provides the cell-cell communi¬ 
cation. X then excites a downstream species while Y 
inhibits it (LEGI). Conceptually, X measures the local 
concentration of external molecules, while Y represents 
their spatially-averaged concentration. If the local con¬ 
centration is higher than the average (i.e., the excitation 
exceeds the inhibition), then the cell is at the higher edge 
of the gradient. While such comparison of the excitation 
and the inhibition can be done by many different molec¬ 
ular mechanisms |13j . here we are interested in the limit 
of shallow gradients. In this limit, biochemical reactions 
doing the comparison can be linearized around the small 
difference of X and Y, and the comparison is equivalent 
to subtracting Y from X [4j. Therefore, we take this 
difference, A, as the readout of the model. 

Because we are interested in the limits to sensory pre¬ 
cision, we focus on the most sensitive regime, the linear 
response regime, where the effects of saturation are ne¬ 
glected. Introducing r„, and 2 /„ as the molecule num¬ 
bers of R, X, and Y in the nth cell, the stochastic model 
dynamics are 

N 

C = DV'^C- ^ 6{X-Xn)fn, 

n—1 

rn = aCn - + rin, 

Xn = Prn - VXn + ?n, 

N 

Vn U ^ ^ Xni (3) 

n' — l 

where c{x, t) is the external concentration, and c„ = 
c{xn , t) is the concentration at the location of the nth 
cell. The matrix M„„/ = {l + 2"f/i^)Snn' - ("f/i^)iSn',n-i + 
^n',n+i) includes the neighbor-to-neighbor exchange 
terms and is appropriately modified at the endpoints 
n G {1, A}. r]n, and Xn are the noise terms. Specif¬ 
ically, arises from the equilibrium binding and un¬ 
binding of external molecules to receptors and can be 
expressed in terms of fluctuations in the free energy dif¬ 
ference Fn associated with one molecule unbinding from 
the nth cell [14], ijn = aCnSFn (here Fn is in units of the 
Boltzmann constant times temperature). The Langevin 
terms and Xn account for noise in the activation, deac¬ 
tivation, and exchange reactions. They have zero mean 

and obey m 

= Sn'nWfn + VXn)5{t - t'), 

{Xn{t)Xn'{t')) = [5n'n{Pfn + I’Vn + ‘^IVn + IVn-l + IVn+l) 

- 5n'+ IVn) 

- 5n',n+l{iyn+l+iyn)]5{t-t'), (4) 

where positive terms account for the Poisson noise cor¬ 
responding to each reaction, and negative terms account 
for the anti-correlations introduced by the exchange. We 
are particularly interested in the SNR for the difference 
between local and global molecule numbers in the edge 
cell, Aat = xm — yN, which is the analog of Eqn. [^for 
the idealized detector. 
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A is given by the means Xn and um, which follow from 
Eqn. in steady state: Xn = {j3/v)fN = [a^/(/«^)]cAr 
and vn = {P/v)Y.u^Nlj'n = [a/3/(Mf^)]E„^WnCn, 
such that 


Aat = — 



^n^N—n 


(5) 


Here is the communication “kernel”, 

which determines how neighboring cells’ concentration 
measurements are weighed in producing the global 
molecule number in the edge cell. Previously we showed 
[3] that Kn is 
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To find the noise, we calculate the power spectra of 
xn and j/AT. As explained below, we assume that the 
measurement integration time T is longer than the re¬ 
ceptor equilibration time (ti), the messenger turnover 
time (r 2 ), and the messenger exchange time (ra). Un¬ 
der this assumption, covariances in long-time averages 
are given by the low-frequency limits of the power spec¬ 
tra, = lim^^o[-S'“^/(w) = {5x*Jy^,)]lT. Linearizing 
Eqn. ^ around its means and Fourier transforming (de¬ 
noted by ~) in time and space obtains 

—iuiSc = —Dk^Sc + iuj ^ , 

n 

— iujSXn = aSc{Xn,Uj) — ySrn + aCnSFn, 

-ilxSXn = /3Srn - vSXn + in, 

-iuj6yn= PSrn- ly'^Mnn'SVn’ +Xn, (7) 

7l' 


where 5c{x,oj) = f d^k (27r)~^Sc(k,oj)e~^^'^. The first 
step in finding the noise is to calculate the power spec¬ 
trum for r„, which we do using the fluctuation-dissipation 
theorem (FDT) as in [T3]. FDT relates the power spec¬ 
trum S'”” / (w) (fluctuations) to the imaginary part of the 
generalized susceptibility G'„„'(a;) (dissipation), 

5'w(w) =-Im[Gw(w)], ( 8 ) 

OJ 

where Gnn' (w) describes how the receptor binding relaxes 
to small changes in the free energy. 


Svn = ^ Gnn'{io)5~Fn' ■ (9) 

n' 


We solve for Gnn'ioj) by eliminating Sc from the first 
two lines of Eqn. which yields a relationship between 
Svn and (5F„. As detailed in Appendix [a| writing this 
relationship in the form of Eqn. requires inverting a 
Toeplitz marix (a matrix with constant diagonals), which 
has a known inversion algorithm [20) . The result is 


5 ;;. (a;) 


‘^CtCn' f + ^TTao) n — n, 

1 n 1 — Et n' ^ n. 

^ y A-KaD \n—n \ ' 


( 10 ) 


Here the cell diameter a appears because we cut off the 
wavevector integrals at the maximal value k ~ tt/o, as 
in m- This regularizes unphysical divergences caused 
by the d-correlated noises in the Langevin approxima¬ 
tion in Eqn. In deriving Eqn. |10[ we have made the 
first of our timescale assumptions, namely T ^ ri = 
yT^ + K/FnaD, where K = a/y is the equilibrium con¬ 
stant, and a = a/2 is the cell radius. ri is the receptor 
equilibration timescale: it is the time it takes for a signal 
molecule to unbind from the receptors and diffuse away 
from the cell into the bulk m- Its first term is the intrin¬ 
sic receptor unbinding time, and its second term accounts 
for rebinding events before the molecule diffuses far away 


The second step is to calculate power spectra for ccjv 
and yj\r using the last two lines of Eqn. 


QXX 

^NN 


(w) = 


SIW = MnuM-I, [/?25L-(a;) + {Z^Xw)] 


^NNi^) — 


nn' 

1 


v{v iui) 


( 11 ) 


where M„„/ = Mnn' — iiy:lv)Snn' ■ Now taking the low- 
frequency limit imposes our second timescale assumption, 
namely T 3> T 2 = 1/iz, where T 2 is the timescale of mes¬ 
senger turnover by degradation. The noise spectra in 
Eqn. 11 follow directly from Fourier transforming Eqn. 
and using the steady state means of Eqn. [^to eliminate 
Zn, 


= ‘2i'Xm, 
iXnXn'} ~ 


( 12 ) 


The appearance of M„„/ in Eqn. is expected, since 
the noise arises in reactions in every cell, and then prop¬ 
agates to other cells via the same matrix as the means 
[23) . Indeed, this simplifies the second line in Eqn. 11 
for w —>■ 0, since then M„„/ = M„„/(u; = 0) hits its 
own inverse. The result is an expression for the vari¬ 
ance {SAn)^ = (Sxn)'^ + {SyN^ - + 

‘S'ArA(O) - 2S'^^^(0)]/T, namely 




- 2Y,Kn-, 


SWNiO) , 

T 

SZiO) 


kn-uKn-t 


T 


T 


+ {xn + Koy^). (13) 


Eqns. and [T^ together with Eqns. and [T^ give the 
SNR = (Aat/^Aat)^, which we do not write here for 
brevity. 

The SNR is compared with the result for the ideal¬ 
ized detector (Eqn. in Fig. We see that whereas 
the SNR for the idealized detector increases indefinitely 
with the number of cells N, the SNR for the model with 
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FIG. 2: Precision of gradient sensing with temporal 
integration. Signal-to-noise ratio (SNR) vs. number of cells 
N is shown for the idealized detector (Eqn. with prefactor 
I/tt ) an d for our model with communication (Eqs. [6 10 
and |13[ ). Whereas the SNR for the idealized detector increases 
indefinitely, the SNR for the model with communication satu¬ 
rates for N ^ no. The saturation level is bounded from above 
by the fundamental limit, Eqn. As shown, the bound is 
reached in the high-gain regime a/a^n = fijv = 100 , where 
intrinsic noise is negligible. Other parameters are a = 10 /rm, 
cjv = 1 nM, g = 1 nM/mm, D = 50 fj, = ty = 1 s~^, 

no = = 10 , and the integration time scale is T = 10 s. 


communication and temporal integration saturates, as in 
the no-integration case |4] . This is our first main finding: 
communication leads to a maximum precision of gradi¬ 
ent sensing, which a multicellular system cannot surpass 
no matter how large it grows or how long it integrates. 
The reason is that communication is not infinitely precise 
over large lengthscales. In the next section, we make this 
point clear by deriving a simple fundamental expression 
for the maximum value of the SNR. 


C. Fundamental limit to sensory precision 


The saturating value of the SNR is obtained in the 
limit of large N. In this limit, and when communica¬ 
tion is strong (7 >> v), the kernel (Eqn. re duces to 
K{n,no) « jno [4]. Here uq = y'^jv sets the 

lengthscale of the kernel and therefore sets the number 
of neighboring cells with which the edge cell effectively 
communicates. The limit 7 3> and our assumption 
T = X/v imply our third timescale assumption, 

T T-j, = 1/7, i.e. that the integration time is longer 
than the timescale of messenger exchange from cell to 
cell. Inserting the expression for K{n,no) into Eqn. 
and approximating the sum as an integral in the large N 
limit, the mean becomes A^y ~ (pcfd/— Cm-uq) = 
afinoag /fjiu. Inserting K{n,no) into Eqn. 13 results in 


products of the exponential with the l/|n — n'\ depen¬ 
dence of the bound receptor power spectrum (Eqn. 10 1 , 


leading to sums of the type e ^/j, which we evaluate 
Appendix |B| The result is 


1 




N 


SNR 

Ceff = Cat -b 


> 




Aat J Tr(noag)^aDT 

log no 


, where 


2 nn 


-(cAr-no/2 “ 2 ca/), 


(14) 

(15) 


and T ^ {ri, T 2 , ra}. Eqn. 14 is fundamental in the sense 


that it does not depend on the details of the internal sen¬ 
sory mechanism. Rather, it only depends on the proper¬ 
ties of the external signal {c,g,D), the physical dimen¬ 
sions (a), and the fact that information is integrated (T) 
and communicated (no) by the cells. The inequality re¬ 
flects the fact that the righthand side contains additional 
positive terms arising from the finite number of bound 
receptors and intracellular molecules (see Appendix [B| . 
These terms represent intrinsic noise and can in principle 
be made arbitrarily small by increasing the gain factors 
a/a^g and jd jv, which dictate the internal molecule num¬ 
bers. What remains in Eqn. |14|is the communicated ex¬ 
trinsic noise, which arises unavoidably from the diffusive 
fluctuations in the numbers of the ligand molecules being 
detected. Eqn. [T^ is shown to bound the exact SNR in 
Fig. HI 

Comparing Eqn. |T^ to the expression for the ideal¬ 
ized detector (Eqn. 0, we see that the expressions are 
very similar but contain two important differences. First, 
whereas Eqn. [2] decreases indefinitely with N, Eqn. [^re¬ 
mains bounded by no for large N (see Fig. [^. Evidently, 
a very large detector is limited in its precision to that of a 
smaller detector with effective size uq. This limitation re¬ 
flects the fact that reliable communication is restricted to 
a finite lengthscale. Importantly, Eqn. demonstrates 
that this noise is present independently of the number of 
intrinsic signaling molecules in the communication chan¬ 
nel. Thus the fundamental sensory limit is affected not 
only by the measurement process (as in BP theory), but 
also unavoidably by the communication process. 

The second important difference is that Eqn.[2]depends 
on Cat, whereas Eqn. 14 depends on the effective concen¬ 


tration Ceff, defined in Eqn. |15| CeS is a sum of the con¬ 
centration measured by the local species, by the global 
species, and the covariance between them, respectively. 
The local species measures the concentration only within 
its local vicinity, cm- However, the global species effec¬ 
tively measures the concentration in the vicinity of ng 
cells. This fact reduces the noise associated with this 
term (and the covariance term) by the factor of ng in the 
denominator of Eqn. |15[ Because intercellular molecular 
exchange also competes with extracellular molecular dif¬ 
fusion, not all of the measurements made by these ng cells 
are independent. Therefore, the reduction is tempered by 
the log ng factor in the numerator of Eqn. [T^ This log 
arises from the interaction of the e“” exchange kernel 
with the l/|n — n'| diffusion kernel (Appendix . The 
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net result is that, because of the correlations imposed by 
externals diffusion, the number of independent measure¬ 
ments grows sublinearly with the system size. Nonethe¬ 
less, it does grow, and, correspondingly, in Eqn. the 
measurement noise in the global species decreases with 
the communication lengthscale ng. Crucially, this means 
that Eqn. is dominated by the measurement noise of 
the local species, i.e. the first term in Eqn. 

In deriving the precision of gradient sensing, we have 
also derived the precision of concentration sensing by 
communicating cells. Specifically, by focusing only on 
the global species terms, and following the steps leading 
to Eqn. we get 


where Oeff = ang/logno and Ceff = c^_„p/cAr_„^/ 2 . 
This expression has the same form as the BP limit, 
(i5c/c)^ ^ l/(aDTc). Indeed, in the absence of a gra¬ 
dient, Cn = Cn is constant, and CeS —>■ cn- Importantly, 
however, the effect of communication remains present in 
Ueff: messenger exchange expands the effective detection 
lengthscale by a factor rig, while ligand diffusion once 
again tempers the expansion by logug. The net result 
is that communication reduces error by increasing the 
effective detector size, Oefi > a. 


SyN 

Vn 


> 


1 


27raeffT>Tceff ’ 


(16) 



FIG. 3: Regional excitation—global inhibition (REGI). 

Signal-to-noise ratio (SNR) is enhanced by allowing both mes¬ 
sengers to be exchanged between cells. The optimal enhance¬ 
ment over LEGI is substantial and occurs because exchange of 
the local species reduces measurement noise, but also reduces 
the signal. Parameters are as in Fig. but with N = 100, 
a/a^y, = P/v = 5 and several values of Uy = \fyyjv as indi¬ 
cated. 


D. Optimal sensing strategy 


In the previous section, we saw that the limit to the 
precision of multicellular gradient sensing is dominated 
by the measurement noise of the local species in the edge 
cell (Eqns. 14 and 15). In contrast, the measurement 
noise of the global species is reduced by the intercellular 
communication. This finding raises an interesting ques¬ 
tion: could the total noise be further reduced if the local 
species were also exchanged between cells? To explore 
this possibility, we extend the model in Eqns. and 
to allow for exchange of the local species at rate 7 ^,, and 
we take 7 —t 7 y > 7 a; for the global species. An imme¬ 
diate consequence of this modification is that the signal 
becomes Ajv ~ 0 iP{ny — nx)ag/yv, where rix = ‘Ixjv and 
Uy = 7 j,/iz. Thus the signal is reduced by local exchange, 
since increasing 73 , decreases the difference Uy — Ux- This 
is because the signal is defined by the difference between 
the global and local readouts, and allowing for the local 
species exchange makes the two readouts less different. 
We thus anticipate that any useful local exchange rate 
will satisfy ^x ^ ly to maintain sufficiently high signal. 
In this limit, we find that Eqn. remains dominated 
by the first term, even as local exchange reduces this 
term according to Cn —t CNi}ognx)/‘^Ux- Eqn. 14 then 
becomes 


For large n^, this expression has a maximum as a func¬ 
tion of nx- The maximum arises due to a fundamental 
tradeoff: exchange of the local species reduces the sig¬ 
nal, but it also reduces the dominant local measurement 
noise. 

The optimal value n* depends on Uy. Experiments in 
epithelial cells suggest that the communication length- 
scale Uy is on the order of a few to ten cells [1]. In this 
range, we find the optimum numerically from the exact 
SNR, which comes from straightforwardly generalizing 
Eqn. [T^ (Appendix [C|). Fig.j^shows that u* is about half 
of Uy for one specific set of parameters, leading to an op¬ 
timal exchange rate ratio of "i^hy = {n%lnyY ~ 25%. 
While the exact optimal ratio depends on the relative 
strengths of different noises, and hence on the gains (see 
Fig. i in Appendix the main finding is robust: a 
multicellular system should exchange both antagonistic 
messenger molecules, one at a fraction of the rate of the 
other. We call this strategy regional excitation-global in¬ 
hibition (REGI). Fig.|^shows that the enhancement over 
the one-messenger LEGI strategy can be substantial. For 
example, with Uy = 10 cells, the SNR is optimally en¬ 
hanced by a factor of 5. With Uy = 15, the enhancement 
is almost 8 -fold. 


II. DISCUSSION 


SNR < (n^ - nx)^nx 

^ CN log Ux 


( 17 ) 


Gellular sensing of spatially inhomogeneous concentra¬ 
tions is a fundamental biological computation, involved 
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in a variety of processes in the development and behavior 
of living systems. Like binocular vision and stereophonic 
sound processing, it is a process where the sensing is 
done by an array of spatially distributed sensors. Thus 
the accuracy of sensing is limited in part by the phys¬ 
ical properties of the biological machinery that brings 
together the many spatially distributed measurements. 
Understanding these limits is a difficult problem. 


Here we solved this problem in the case where the com¬ 
munication is diffusive, and one-dimensional distributed 
measurements are used to calculate external concentra¬ 
tion gradients within the LEGI paradigm. We allowed 
for temporal integration, extending the results of Ref. [1] . 
Some of the features of the gradient sensing limit we de¬ 
rived (Eqn. 14), such as the unbounded increase of the 


SNR with the diffusion coefficient of the ligand or with 
the integration time, carry over from the Berg-Purcell 
theory of gradient sensing [T^], which does not account 
for communication. However, our most important find¬ 
ing is that, in contrast to the BP theory, the growth 
of the sensor array beyond a certain size stops increas¬ 
ing the SNR. The effect is independent of the intrinsic 
noise in the communication system and thus represents a 
truly fundamental limitation of diffusive communication 
for distributed sensing. In particular, it holds for mul¬ 
ticellular systems, as well as for large individual cells. 
Although we derived the limit for a linear signal pro¬ 
file, we anticipate that the limit for a nonlinear profile 
will be similar, Eqn. but with a different effective 
concentration Ceff. It remains to be seen if similar lim¬ 
its hold when the sensors are arranged in two- or three- 
dimensional structures, or when concentration informa¬ 
tion propagates super-diffusively, as is possible in wave- 
based or Turing-type models of polarization establish¬ 
ment [18) . Additionally, it will be important to relax 
various assumptions of the model, such as allowing for 
saturation of receptors or limiting the total number of 
messengers, and coupling the model to the motility ap¬ 
paratus to investigate how the improved sensory precision 
affects downstream functions. 


Our results illustrate two important features of tem¬ 
poral averaging by distributed sensors. First, our deriva¬ 
tion naturally reveals which timescales are relevant in 
this process, namely receptor equilibration (ri = + 

K/ATraD), messenger turnover (t 2 = Ajv), and messen¬ 
ger exchange (ts = I/ 7 ). In principle, these timescales 
could have depended on system-level properties, such 
as the system size {N) or the communication length 
(no). Surprisingly, instead they depend only on single¬ 
cell properties, meaning that efficient temporal averaging 
is not slowed down by increasing the number of sensors. 
Second, our results reveal the effects of over-counting due 
to correlations between external and internal diffusion. 
In both gradient sensing (Eqns. 14 and 15) and concen¬ 
tration sensing (Eqn.[I^, we see that the noise reduction 
afforded by communication-based averaging is tempered 
by a factor log no. This log is not a mathematical cu¬ 
riosity. Rather, it reflects the fact that not all measure¬ 


ments communicated to a cell by its neighbors are inde¬ 
pendent since the signal molecules also diffuse externally. 
Coupling external diffusion with internal exchange intro¬ 
duces correlations among measurements, which reduces 
the benefit of internal averaging. 

Another central prediction is that the gradient sens¬ 
ing is improved by a system with two messengers, ex¬ 
changed at different rates. We call this mechanism re¬ 
gional excitation-global inhibition (REGI), a generaliza¬ 
tion of the standard local excitation-global inhibition 
(LEGI) model. Optimality of REGI follows directly from 
the interplay between the ligand stochasticity and the 
communication constraints. Therefore, REGI has not 
been identified as an optimal strategy in previous stud¬ 
ies that neglected either of these two effects. However, 
evidence for REGI may already exist in large gradient¬ 
sensing cells, where activated receptor complexes, which 
diffuse in the membrane at a rate ^ 10—100 times slower 
than similar cytosolic molecules [H], may act as the re¬ 
gional messengers. 

REGI emerged from maximizing the SNR in our sys¬ 
tem, which revealed the optimal rate ratio '^x/ly Maxi¬ 
mization of the SNR also implies that the optimal value 
of 7 y (or ny) is infinity, since the SNR grows indefinitely 
with Uy (see Fig. Infinite riy corresponds to averaging 
over as large a distance as possible. Such a strategy is 
only optimal here because the concentration profile is lin¬ 
ear, with constant gradient g. In contrast, more physical 
nonlinear profiles (e.g. exponential, power-law, randomly 
varying, or profiles with extrema) have spatially vary¬ 
ing g. In these cases, if the size of the group of cells is 
larger than the correlation length of g, then an infinite 
Uy would average out the signal together with the noise, 
which would reduce the SNR. In contrast, a finite Uy 
would allow a subset of cells to detect the local gradient 
in their vicinity, which is an essential task in morpho¬ 
logical processes such as tissue branching and collective 
migration. 

REGI can be interpreted as performing a spatial 
derivative. Specifically, the two-lobed filter K{n,nx) — 
K{n,ny) reports the difference in concentrations mea¬ 
sured over distances rix and Uy near a given detector, 
and the values Ux and ny depend on the properties of 
the environment. Thus REGI is similar to the temporal 
differentiation in E. coli chemotaxis. Indeed, the tempo¬ 
ral filter of the E. coli sensory module is also two-lobed, 
with the short and long timescales set by ligand statis¬ 
tics and rotational diffusion, respectively [53]. Thus, for 
both spatial and temporal filtering, the choice of the two 
optimal length- or timescales is determined by matching 
the filter to the statistical properties of the signal and 
the noise [2^, which is understood well for E. coli [25] . 

Interpreting the REGI model as a spatiotemporal fil¬ 
ter suggests experiments that would identify if a certain 
biological system employs this mechanism. Such exper¬ 
iments would involve concentration profiles that differ 
substantially from steady-state linear gradients. For ex¬ 
ample, subjecting cells to a concentration profile with a 






spatially localized maximum would allow one to measure 
both 7 a; and jy by observing the response of cells near 
the concentration peak as a function of the peak width. 
Alternatively, one can subject cells to a spatiotemporally 
localized concentration pulse and observe if a response a 
certain distance away from the pulse exhibits the signs of 
only inhibition (LEGI, one messenger), or inhibition and 
excitation on different scales (REGI, two messengers). 
Understanding fundamental sensory limits for diffusive 
communication in gradient sensing opens up possibilities 
to propose and analyze these and other related experi¬ 
ments. 
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Appendix A: Power spectrum of the bound receptor number 

The dynamics of receptor binding are given in Fourier space by the first two lines of Eqn. 

—iuiSc = —Dk^Sc + iuj ^ , 

n 

— iuiSrn = aSc(Xn, cj) — flSrn + aCnSFn, 

where 

Sc{x,ijj) = [ d^k{2TT)~^5c{k,oj)e~''^'^. 


We solve Eqn. A1 for 5c and, using Eqn. A3 insert it into Eqn. 


A2 


to obtain 


- iw[l + aY.{uj)]}5rn - iuia ^ F(|f„ - f„/|,w)(5r„/ = aCnS'Fn 

n'^n 


(Al) 

(A2) 


(A3) 


(A4) 


where 


S(cc) = 
V{x,uj) = 


(Pk 


(27r)3 Dk'^ — iu) 

d^k _ 

(27r)3 Dk'^ — ioj 2 'k‘^x 


1 

27r^ 

1 


dk 


dk 


Dk"^ — iuj ’ 
k siii{kx) 


and 


Dk'^ — iu 


(AS) 

(A6) 


are the “self-ener gy” and “int eraction potential” between cells mediated by diffusion, respectively [Hj. The simpli¬ 
fications in Eqns. AS and A6 come from writing the volume element d^k = k^ sin 0 dk d0 d(j) in spherical coordinates 
aligned with x, such that k ■ x = kx cos 0. 

We are interested in the low-frequency limits of E(a;) and V{x,uj). V{x,0) = l/(47rl?x) is finite, whereas S(0) 
diverges. The divergence stems from the delta functions in the dynamical equations, which model the cells as point 
sources. As in [14], we regularize the divergence by introducing a cutoff A ^ tt/o at large k to account for the fact 
that cells have finite extent a, making E(0) ~ (27r^)“^ dk/D = l/(2T:aD). This models cells as spheres of diameter 
a, but the exact shape of the cell will not be important for the limits we take. These expressions allow us to write 
Eqn. jArjas 


N 

^ ( ^nn 


■Sr„> = 


aCri 




(1 -I- iujTi) S~Fr. 


where 


Tl = 


1 a/u 
-+ ' 


fj, 2'KaD ’ 


= 5n 


Z = -lUJ- 


5nn0 




47raZl 


(A7) 

(AS) 

(A9) 

(AlO) 


In writing Eqn. A7 we have assumed that ujti is small. This assumption is valid for integration times T = 27r/a; 
much longer than ri. The quantity ri is the receptor equilibration time: it is the time it takes for a signal molecule 
to unbind from the receptors and diffuse away from the cell into the bulk [21j . Its first term is the intrinsic 
receptor unbinding time, and its second term K/FkgD (where K = a/ is the equilibrium constant and cr = a/2 is 
the cell radius) accounts for rebinding events that occur before the molecule diffuses away from the cell completely 
[22| . Either term can dominate: the first term dominates if the intrinsic association rate a is much smaller than 
the diffusion-limited association rate FkgD, because then the molecule rarely rebinds. Conversely, the second term 
dominates if a is much larger than 47rcrll, because then rebinding is frequent, and comprises most of the escape time. 
We see from E qns. A8 and AlO that \z\ < ujti; therefore, we also treat 2 as a small parameter. 

Solving Eqn. A7 for dr„ requires inverting the matrix This matrix is a Toeplitz matrix (a matrix with constant 

diagonals), which has a known inversion algorithm |20] . Since Ann' is also symmetric, it is completely specified by its 
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first row \\ pi p 2 ... pn-i]- The inversion is performed recursively as follows. First one introduces N — \ scalars hi, 
/i 2 , ..., hN-i and N — 1 column vectors ..., These are initialized as 

hi = l-{pif, g<') = [-pi], (All) 


and updated as 




hk 


^k+i) 


Ck/hk 


(A12) 


where (^k = Pk+i + qisq with the elements in reverse order. Then the inverse is written in terms of 

the final quantities /itv-i = H and = Q. The upper left element is 

An = ^/H, (A13) 


the rest of the first row and column are 


Ai+1,1 — — Qi/H (1 < * < 1), 

and the diagonals are calculated recursively from the first row and column as 

-1 _ r-i 


(A14) 


L-^i,,+i=L-/ + {Q.Q,-Q,Q,)/H {1 < {i,j} < N - 1). (A15) 

We now apply this algorithm to Eqn. |A9[ keeping only terms up to first order in the small quantity z. From Eqn. 


A9 we have pj = zjj. The initial values are 

/ii = 1 — « 1, 

Using ^1 = z/2 + (z)(z) « zj^, the first recursive step gives 


qTD = [-z]. 


/l2 = 1 - 


1 




r-z-(z/2)(-z)/(l)] 


—z 

[ -(^/2)/(l) J 


-z/2 


Continuing the recursion establishes the general formulas 


hk = 1 , 


AC 


from which we extract the final values, 


AT = 1, 


= -zjj (j < k), 


Q] = -^/j- 


(A16) 

(A17) 

(A18) 

(A19) 


These values immediately provide the first row and column of the inverse according to Eqns. |A13| and |A14| Then 
noting that the second term on the right hand side of Eqn. |A15| is second order in z, that equation implies that the 
diagonals of the inverse are constant. We therefore have the inverse to first order in z, 


A„ri' “ ^rin' 


^nn') 


\n — n' 


The inverse allows us to solve Eqn. A7 for (5r„, 

Srn = ^ Gnn'5'Fr, 


N 


(A20) 


(A21) 


where 


Gnn'{<^) = - (1 + = 


a.Cn 




ATTaDfi'^ \n—n' \ 


2'KaD) 


n = n, 


is the generalized susceptibility. The fluctuation-dissipation theorem then gives the power spectrum, 


0 ^ eye 

SZA^) = -Im[Gw(a;)] = 

CJ 


'(1 


2'KaD) 
1 


4:KaD \n—n'\ 


n = 71 , 

v! yf n, 


(A22) 


(A23) 


as in Eqn. 
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Appendix B: SNR in the many-cell, strong-communication limit 

The variance of the readout is given by Eqn. [T^ 


{SAn)"^ = 




■’ NNy 

T 


n^n' — l 


'' Nn\ 

T 


+ [xn + KoUn) , (Bl) 


where xn = af3cN/tJ'i^ and ijn = Q-P X]ti=o KnCN-n/IJ- v. In the limit of many cells {N ^ 1) and strong communication 
(7 ^ v), the kernel takes the approximate form 


K ^ p-nlnp 
no 

where ng = \/^jv S> 1 is the communication lengthscale. Using Eqns. 
term. 

The first, fourth, and fifth terms in Eqn. |Bl]are straightforward, 


(B2) 


A23 


and 


B2 


^2,5-^(0) 2 a/32 1 ^2^2 _ 


1/2 T 


/rT /rz/2 ' naDT 

2 a/3_ 

tv;—C i\ 

'1 

2 a/3 cjv—no 


Cat, 


= - 7 ^ -Cat, 

vl vl fiiy 


— KoVn = „ 

vl vl fiv Uq 


we evaluate Eqn. B1 term by 

(B3) 

(B4) 

(B5) 


In the last equation we use the limit of large N to approximate the sum in ypf as an integral, 

N-l 


i rOO 2 

KnCN-n « / dn —e“"/"“(cAr - nag) = cn - 

t'o do rio 


^0^9 — ^N—riQ' 


The third term in Eqn. |B1| we split in two, 

/^9 N 




n—1 


T 1/2 


Kc 


5]vV(0) 


N-l 




n—1 


' T 


The first of these is like Eqn. |B3| 


_2^K ^]vV( 0 ) _ 4 a /32 _ 2 a2/32 c^y 

z/2 ° T gT /rz/2 no naDT gPv^ no 


The second evaluates to 

2 N-l 


n—1 


2 N-l 

W r-iN-n)/no 

T z/2 '' AiraD N — n 

n—1 


- 2^ = -2g 4|.[c-„ - (iV - n)a„l “ ^ 


1 a2/32 

naDT gPv'^ 

1 a2/32 


- Af-l _,/nn ^-1 

^ e _<Pl'y g-j/no 

”0 •?' ^ 


Cjv , 

— log no - ag 
no 


naDT fjPv^ 

where the last step again uses the integral approximation for large N. In particular, we have written 


N-l 

E 

3=1 


.-3/"o g-3/"o 


dj 


= r(0, l/no) « log no - 7 e « log no. 


(B 6 ) 


(B7) 

(B 8 ) 

(B9) 

(BIO) 

(Bll) 


(B12) 
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using the small-argument limit of the upper incomplete Gamma function. The last step assumes that log no is much 
larger than the Euler-Mascheroni constant 7 e « 0.577. Numerically, we have checked that Eqn. |B12| is valid for 
5 ^ Altogether, the third term in Eqn. B1 is then 




■’ Nny 

T 


4 a/3'^ cj\[ 2 a^/3^ 

/in^ no naDT [no 2no 

4 Cm 2 a^/3^1ogno 

/iT no iraDT 2no 


Cat Cat ag 

-^ Tir - y 


-Cat, 


(B13) 

(B14) 


where the second step uses cm — agno/2 = CM-no /2 < cm and assumes log no 2. 
The second term in Eqn. |B1| we also split in two. 


^ ^N-n^N-n' -= ^ 


err err /n\ 

A^-n- If -f 2^ 2^ KM-n^N-n' -- 


n—1 


n—1 n'^n 


The first of these is straightforward to evaluate with the integral approximation in Eqn. |B6 


13'' 


,,2 Z^^N-r 


S^niO) _ 1 ^_,^/^„2acM- 


n—1 


T 


— tL. __ 

-2^ qri2 


j=0 






2 a/3^ CM-na/2 1 Oi^(3'^ CM-no/2 


fiT 2no iraDT 2no 


The second can be evaluated in two parts 
52 N 
V 


TZ ‘5'nn'(0) 
l3M-nKM-n' - Tf, - 


1 


n—1 n'^n 


2T:aDT 


_ N-l 

Cm 


g-{j+f)/no f^g j'g“(t+t')/"o 


CM gw.-,,- oiy 

77 2 I 7 _ 7/1 ^2 2-^ 


ns ^ —' 17 — j’ 

° 7=0 7 V 7 


ns —' I? — 


(B15) 

(B16) 

(B17) 

(B18) 


A B 

Notice that -B = — i9„A/2, where u = 1/no, so that we only need to evaluate A. We split A into two equal components. 


22 22 e~^3+j')/n.o 22 22 e~^^~^3')/'no 

^ r- + i: E 


j - j 


22 g-(7+7')/rao 


7=0 7'=0 
TV-l 7- 

= 2 EE 


7=0 7'=7 + l 


J - J 


7=0 7'=0 

and rewrite it in terms of fc = j -I- j' and i = j — j' 


J-J 


(B19) 

(B20) 


AT-i 2(7V-r)-f _J,, 

^ = 2E E V' 

l=\ k=i,i+2,i+i,... 


(B21) 


We approximate with integrals for large TV, accounting for the fact that k has support on only half of the integers in 
its range. 


rN 


A 


di 




dk 


= no d£ 


N ^-i/no 


- noc 


-2N 


rN pt./no 

d£ 2_. 


(B22) 

(B23) 


The first integral is approximately log no by Eqn. B12 The second integral evaluates to Ei(A^/no) — Ei(l/no), where 
Ei is the exponential integral function, whose large- and small-argument limits are Ei(A^/no) « /[N/no) and 

Ei(l/no) ~ — log no, respectively. Thus, the second term in Eqn. B23 vanishes exponentially with iV, and we have 


A = no log no, 

^ = y(l+ logno) « y logno. 


(B24) 

(B25) 
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second term in Eqn. B1 is then 


making the term in brackets in Eqn. B18 equal to {cn — agno/2){logno)/no = logno)/no- Altogether, the 


N 


^2 ^N-uKn-t 

n^n' — l 


T 


2 aj5^ CN-no/2 
/iT 2no 


1 a^/3^ 

iraDT 


CN-no /2 log no 
2 no 2 no 


CN-no/2 


2 aP'^ CN-no/2 1 a^/3^1ogno_ 


/lT 2 no 


naDT 2no 


CAr-no/2i 


(B26) 

(B27) 


where the second step assumes log no » 1 . 

Einally, collecting the terms in Eqns. |B3||B5[ |B14[ and|B27[ the variance in Eqn. |B1| becomes 




a/3 

a^/jiv 


- , logno_ 1 3 

ttOT “ 


Extrinsic noise from c, propagated to x and y 

2 a^/x f _ 
fjT a 


Cn 


CN-no/2 _ 2— ) 


2 no 


no 


Intrinsic noise in r, propagated to x and y 

2 n 3 2 v (CN-no \ 3 

vl a p vl a p \ no 


(B28) 


Intrinsic noise in x 


Intrinsic 


The last line contains the intrinsic noise from xn and ijN- For example, the first term in this line is = 

{2/vT){apCN /= 2.xn/'>^T] the relative noise [5xn/xnY — ‘^■/vTxn then decreases with the number of molecules 
Xn X vT that are turned over in time T, as expected for intrinsic counting noise. The second term in this line is for y 
and is similar, except that since the global species is exchanged over roughly uq cells, more molecules are counted and 
the noise is reduced by a factor uq. The second line contains the intrinsic noise in r, propagated to x and y. The third 
line contains the noise in c, propagated to x and y, which we deem extrinsic, since it originates in the environment 
and is not under direct control of the cells. 


Importantly, the intrinsic noise terms in Eqn. B28 are reducible by increasing the numbers of receptors and local and 
global species molecules. These molecule numbers are set by the gain factors tnIcl^cn = a/a^p and xn/^n = /3/x^j 
and indeed, we see that the second and third lines in Eqn. |B28| vanish as the gain factors grow large. In this limit we 
are left with only the lower bound set by the extrinsic noise, 


\2 ^ f 1 f- , log? 


(B29) 


Dividing by the square of the mean An = ajSnoag/pv produces Eqns. 14 and 15 


Appendix C: Exact SNR for regional excitation—global inhibition (REGI) 


In the REGI strategy, both messengers X and Y are exchanged between cells, at rates px and 7 ^, respectively. This 
results in a straightforward generalization of the expression for the signal-to-noise ratio (SNR). Specifically, the mean 
becomes (compare to Eqns. and in the main text) 


A 


N 




"^CN-r, 


N-1 

E 

n—O 


K^CN-r 


(Cl) 


where now there are two communication kernels. 
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FIG. 4: Dependence of the optimal exchange rate ratio on system parameters in the regional excitation- 
global inhibition (REGI) strategy. The signal-to-noise ratio (SNR) has a maximum at a particular rate ratio which 

increases as a function of either gain factor, (A) or (B) /3/i/. Parameters are a = 10 /rm, cjv = 1 nM, g = 1 nM/mm, 

D = 50 /im^/s, T=10s, g = u = l s”'^, N = 100, and Uy = y/'y/v — 10. In A, fijv — 5, and aja?' g, is varied as indicated. In 
B, a/c?g = 5, and Pjv is varied as indicated. 


The variance becomes (compare to Eqn. 131 






—OR- 


N—n^N—n' 


SZu'iO) 

T 


+ — {K^XN+Kym): 


(C4) 


,{u}) remains the same as in Eqn. 10 (or equivalently Eqn. A23). The 


where the bound receptor power spectrum 
SNR is then (Ajv/SAjv)^. 

The SNR has a maximum as a function of the rate ratio ^xhy The location of the maximum ')xllv must lie between 
0 and 1. At ^x/ly = 0, the X messenger is not exchanged, and we recover the SNR of the local excitation-global 
inhibition (LEGI) strategy, which is a limiting case of REGI. At ^xHy = 1, there is no difference between X and Y, 
and the signal (and therefore the SNR) is 0. The exact location of 'fxHy depends on factors that are specific to the 
system, e.g. the concentration profile c„, the environmental and system parameters, and the measurement location 
(n = 1 vs. n = iV). 


We illustrate the dependence of ')xlly on particular system parameters, namely the gain factors a/a^g and P/v. 
For this. Fig. Hshows the dependence on SNR on ^xHy as we vary the gain factors a/a^g (A) and (3/v (B). In both 
cases, we see that the optimal rate ratio ^xHy increases with increasing gain. This is because increasing either gain 
factor increases the number of internal messenger molecules. With more molecules, the system can afford to increase 
7 a; while maintaining the same difference in molecule number Ajv in the iVth cell. The increase in ya, enhances the 
spatial averaging by the X messenger, and thus reduces the noise in the estimate of cat. Therefore, we see in Fig. 
that the maximal SNR occurs at a higher 7 * value as either gain is increased. 



























